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Abstract 

After a derivation of the quantum Bayes theorem, and a discussion of the reconstruction of the 
unknown state of identical spin systems by repeated measurements, the main part of this paper 
treats the problem of determining the unknown phase difference of two coherent sources by photon 
measurements. While the approach of this paper is based on computing correlations of actual 
measurements (photon detections), it is possible to derive indirectly a probability distribution for 
the phase difference. In this approach, the quantum phase is not an observable, but a parameter 
of an unknown quantum state. Photon measurements determine a probability distribution for the 
phase difference. The approach used in this paper takes into account both photon statistics and 
the finite efficiency of the detectors. 
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I. INTRODUCTION 



The use of Bayes theorem in quantum mechanics is discussed. It is shown that the 
quantum Bayes theorem follows from the ordinary quantum measurement theory, when 
applied to density operators that represent our a-priori knowledge of a system. The examples 
studied involve measurements on multiple copies of the same (unknown) state. The theorem 
is used to determine the unknown state by successive measurements on several of the copies 
of the state. 

The theorem is applied to quantum optics: an idealized information-theoretic description 
of propagating CW laser beams is treated in detail. It is shown how photon detections on 
part of the beams can be used to determine the phase of the rest of the beams. Explicit 
expressions are derived for the conditional probabilities of detecting photons at different 
positions, given the numbers of photons detected at different positions in the past. The 
quantitative predictions could be used, in principle, to test models of the quantum state of 
propagating laser beams. 

I hope that the information-theoretic approach used in this paper will be a useful contri- 
bution to the understanding of the problem of quantum phase. 

II. BAYES THEOREM 

Bayes theorem follows from the definition of conditional probability: 



Bayes theorem allows us to reverse probability relationships. For instance, if the probability 
of the observation B, given the assumption (theory?) A is P{B\A), then the probability of 
A, after having observed B is: 



n(y4) is the a-priori probability of A. The form of Eq. Q is the most useful for analyzing 
classical data. The Bayes Theorem requires some way to provide the a-priori probabilities 
(priors). In some cases the priors can be supplied by symmetry arguments, in other cases, 
repeated applications of Eq. (0) will provide a result only weakly dependent on the original 
priors. 
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In quantum mechanics one has to go back to the form of Eq. ^ : 

where 11 is a density operator containing our a-priori knowledge of the system being mea- 
sured. Ab and are measurement operators corresponding to B and A. 

If Ab and Aa commute with each other (for instance if they act on different subspaces — as 
in all the following examples), then Eq. (jS)) can be put in the form 

PiAlB, - »f . (4) 

Here Ea and Eb are the POVMs corresponding to the operators Ab and Aa- 

Eb = A^B^B, Ea = A\Aa. (5) 

Under the same assumptions that led to Eq. one can derive the conditional probability 
of measuring Ai, A2 ■ ■ ■ An, after having measured 5i, i?2 • ■ • B,n- 

TD( A A A \a n a \ "^t^O^Eb^Eb^ ■ ■ ■ Eb^EaiEa^ ■ ■ ■ Ea„) , . 

P{M,A. ■ . . A„|B, . . . B.„) = ______ . (6) 

Note that the quantum Bayes "theorem" as expressed here in Eq. Q and is not a new 
principle, but a consequence of the ordinary quantum measurement theory. The only special 
assumptions that will be made, will be the choices of the a-priori density operators 11. 

III. SPIN SYSTEMS 

An example of the use of Eq. Q is a system of + M qubits (see Refs. Q and 0). 
Following Ref. we assume the a-priori density operator 

flQbits = / d'xn{x) , (7) 



where the vector x has components x, y and z. The vector operator a has components 
ax, CTy and (the Pauli operators). Eq. (|7j) means that all M + qubits are in the 
same state, represented by the same density operator, which is unknown, except for the 
a-priori probability distribution n(x), which is normalized so that its integral is equal to 
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one. The form of Eq. (|7j) follows from the quantum de Finetti representation theorem. (For 
information about the quantum de Finetti theorem see Refs. y and 

The measurement operator corresponding to measuring M^. times cTx (on different qubits), 
My times dy and times cTz (Mj. + My + M^ = M), with the result = +1 obtained 
times and = —1 obtained M~ times (with M^ + M~ = M^), etc. is 

1 ■ ■ ■ ® ( j + Permutations. (8) 

Similarly, the measurement operator corresponding to measuring A^^, times a^, Ny times CFy 
and Nz times (A^^ + A^^^ + A^^ = A^), with the result = +1 obtained times and 
(Tx = —1 obtained A^~ times (with A^+ + A^~ = A^^.), etc. (also all on different qubits) is 

® + Permutations. (9) 

Inserting Eqs. ((Tj), (jH)) and into Eq. Q, we get a result for the probability of measuring 
ax = +1 times, ax = —I N~ times, etc., after having measured ax = +1 times, 
ax = —1 M~ times, etc. (with each measurement performed on a different qubit): 

P(Ar+, N-,N^, N;,Nt, Ar-|M+, M", M+, M;, M+, M") = 

{K + N~y.{N^ + Ny')\(Nt + N;)\ 
N+\N-\ N+\N-\ N+\N~\ 

J.^ ^iy -^'y "^'Z 

±x \ M++N+ ( ]-x \ M-+N^ 

2 I ^ 2 ' ^ 2 ^ ^ 2 ' 2 ' ^2 



l±x \ M++N+ ^ 2-x ^ M-+N^ ^ l±y ^ M++N+ ^ ^ My+Ny ^ l+z ^ M++N+ ^ 1^ ^| M'+N; 
l+x \ M+ I l-x \ M' ( 1+?/ ^ Af + (^-y\ My ( l+z \ M+ (1~z\M: 



.(10) 



2 ' y 2 ' y 2 ' y 2 ' y 2 ' y 2 ' 
Eq. (jlUj) is an obvious consequence of the discussion of Ref. I wrote it down explicitly 
to illustrate the symmetry between Ms and A^s in the numerator. This symmetry has 
consequences for the naive use of the relation 

Af^ooV 2 y V 2 y Mx-^ + Mx' ^ ' 

up to a normalization factor. Eq. (fTT|) is only applicable in Eq. ()10|) if the Ms are large 
compared to the A^s. So a large, but finite number of measurements will determine the state 
accurately enough only if the number of successive measurements is small compared to the 
number of measurements used to determine the state. 
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A. Exact Solutions 



To illustrate this point, we will give a complete treatment of the case in which the spin 
measurements are all along the same axis. We will give an analytical solution with the prior 
of Ref. P]. In this case, the integral of Eq. (fTUI) is constrained to the surface of the unit 
sphere, with uniform weight. (This is — of course — not the only possible pure-state prior: 
one could use a non-uniform weighting of the unit sphere surface.) The a-priori probability 
of measuring times cr^,. = +1 and M~ times = — 1 is 



M+!M-! ./o 7o 47r V 2 



Substituting x for cos^, we can put Eq. (fT^ in the form 



M+'.M-l J_i 2 

The conditional probability of getting times ax = +1 and N~ times ax = —I — after 
having measured times a^ = +1 and M~ times ax = — 1 — is 

(AT+J-M~\\ dx(l±x\M++N+ (]-x\M-+N- 
P{N^,N^\M^,M-) = V — 7^ • (14) 

^ ^ J-1 2^2' ^ 2 ' 

Using the generating function method, one can show that P{M^ ^ M~) has the value 

This means that all possible outcomes of the Mx = M^ + M~ measurements have the same 
probability. (Here we consider equivalent the measurements that give the same number 
of spin ups and downs, independently from permutations — this explains the combinatorial 
factor in Eq. (fT^ and 

Because the integrals in Eq. ()14p have the same form as the integral of Eq. ()13|) , we now 
have all we need to write down explicitly the result for Eq. (fT^ : 

P(N- N-\M- M-) = (^^ + ^^y-(M^ + M^ + mM- + NmMx+N-)l 

N+\N-\ {M+ + M- + N+ + N- + l)\M+\M-\ ' ^ ' 

Eq. (fTBj) is the exact result, valid for any A'^^, A''", M+ and M~. By the use of the approxi- 
mation 

Uo, = (17) 

M^oo Ml ^ 



we can obtain the limit of Eq. ()16|) for ^ A^^: 

,11. P(Ni,M;lMi,M;) . (^Y' . (18) 



M.^oo ^ N+\N-\ \M+ + M-J \M+ + M- ^ 

This is exactly the result that would follow from using Eq. (fTT|) in Eq. (fT^ . We have been 
able to derive it only in the case in which M+ + M~ ^ N^, N~ . 

As a particular example, the probability of getting ax = +1 times, after having 
observed = +1 Af+ times, is: 

P(N+\M+) = ^^—^^ = " + (19) 

(Because for large Mx~^ the integrals in Eq. (fTn|) are dominated by a small neighborood of 
X = 1, if Mx~^ is large, the result of Eq. is actually valid for any prior n(x) that is not 
zero at X = 1.) Now, if ^ N^, the probability given by Eq. (fTU|) is 1, as expected for a 
state that has been "determined" to be = +1. However, if = A^^, then — no matter 
how large is — the probability given by Eq. (fT^ is only i. This shows the difference of 
state "determination" using Bayes theorem and state "projection" from a "von Neumann" 
measurement. A single projective measurement to the state = +1 will force all successive 
measurements of ax to give the result +1, no matter how many the successive measurements 
are. 



This "weakness" of Bayesian state "determination" is of importance 



mation theory because such methods are invoked (for instance in Ref. jsl), to argue that a 
conventional laser can be used for quantum teleportation with continuous variables, contrary 
to recently made claims (see Ref. ^) that teleportation requires novel, truly coherent light 
sources. 



or quantum infor- 



IV. LASER BEAMS 

In the rest of this paper, Bayes theorem will be applied to a model of propagating laser 
beams. We will discuss the case in which the probability of detecting several photons at 
a time is not negligible. The model for the density operator of two laser beams given in 
Ref. y is 

riLase. = r ^ r ? ® |&e^^^)(te'^1)^'' . (20) 

Jo 27r Jq In ^ 
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The prior of Eq. ()2U|) follows from a simple symmetry argument: since there is no reason 
to prefer any particular phase, choosing an uniform a-priory distribution is very reasonable. 
The output beam from each laser is made up by "packages", all with the same, unknown 
phase. (The form of Eq. (j!2(J|) also follows from the quantum de Finetti theorem — see Refs. 
and ||.) 

A coherent state lae*"^) is defined by the eigenvalue condition 



(21) 



where a is a photon destruction operator. I put hats on all operators: all other symbols are 
real numbers (except i, of course). 

We will measure photons in states obtained by combining the two beams in states pro- 
duced by a 50/50 lossless splitter/recombiner. The destruction operators of the states pro- 
duced by the beam splitter are 



and 



d 



d + b 

d — b 
2 ■ 



(22) 



(23) 



The setup, which is identical to that of Ref. 



is illustrated in Fig. ^ The labels c and d 




FIG. 1: The output beams from the two lasers (a and b) are mixed on a beam splitter and the 
resulting beams (c and d) are measured by two photodetectors. 



will be used to denote both the two beams coming out of the splitter and the corresponding 
detectors. The POVM corresponding to measuring photons at c is 



E 



Nc. 



-TjC' C . 

• ) 



(24) 



where the colons represent normal ordering — annihilation operators to the right, creation 
operators to the left. (Time ordering is not necessary, because we are considering a non- 
interacting field — except for the measurement process.) A similar expression gives the 
POVM corresponding to the detection of A^^^ photons at d: 

Here 77 is a constant, dependent on the quantum efficiency of the detectors. The time 
duration of the detection process (which we assume to be short) is included in the constant 
rj. (I will not give a proof of Eq. and ()25|) — the reader will find a proof, for instance, 
in Chapter 14 of Ref. jsf.) Note that the measurement operators of Eq. ()24j) and ()25|) 
consistently take into account the finite efficiency of the measuring devices. To make the 
algebra simpler, we assume the same detection efficiency at the two locations. 

To provide a normalization, if one measured directly the photons from a — by removing 
the beam splitter of Fig. ^ for instance — one would obtain for the probability of detecting 
Na photons 

m) = ^Vr^"' ' (26) 

with average number r/a^. The average number of photons detected directly from b would 
be Tjb^. 



A. Equal Frequencies 

I will use the Heisenberg picture in these sections. This means the the operators Enci 
etc. depend — in general — on time. 

I will begin the discussion by assuming that the frequencies of the two lasers are identical. 
In this case the time dependence of the operators En^ and En^ can be neglected. Inserting 
Eqs. fl20|) ■ ()24|) and ()25p in Eq. ®, we have the probability of detecting Nc photons at c in 
the second package, if photons have been detected at c and photons detected at d 
in the first package: 

P(A,|M„ M,) = Tr(nLase.s(i?M.i?Mj^i?A.J ^ ^2^^ 

Tr ( IlLasers Em^ Em^ ) 



which gives the result 



P{N,\M,,Mi) = ^^^-^ -^^^ — -— ^ . 28 
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In Eq. (j28|) and in the following, denotes the difference of the two phases: cp = (pa — 4>b- 
The integral over the sum of the phases does not affect any of the results, and it will be 
simply omitted. (The results of Eq. and successive equations follow from the fact that — 
because of the normal ordering and the form (j2(J|) of the density operator — one can substitute 
the value <^<^"^°-+^^"*''' for the c operator and the value for the d operator inside the 

integral.) The conditions of Eq. © are satisfied, because the measurements on different 
packages commute, and the measurements at c and d commute, because [d\ c] = [c\ d\ = 0. 
Eq. can be rewritten as 



P{N,\M,,M,) = j^l^ rf0P(0|M„Mrf)(^— + a6cos0)^=e-''("^+'^^™^^), (29) 



where P{(p\Mc, Ma) is the conditional probability of the phase (difference) having the value 
(f). It is given by the expression 

, , , (^ + a6cos0)^^(^-a6cos(/))^^'* , , 



/o d(P'{s^^ + ab cos (PT'^ -ab cos > 
In a similar way, we can derive the probability of detecting A^^ photons at c and 
photons at d in the second package, if Mc photons have been detected at c and Mj, photons 
detected at d in the first package: 

^N. ^N, J^- ^0(a^ + a6cOS0)^^+^^(2^-a6cOS <p)M,+N,^-r,{a^+b^) 



P{N,,Nd\M,,M, 



d 



Nc\Nd\ j^^ dct){^ + ab cos (PY'^i^-ab cos (PY'-i 

(31) 

It is reassuring that one can obtain Eq. ()28|) by summing Eq. ()3ip over A^^. 



1. Low Intensity 

Eq. (j3T| can be specialized to very low photon detection probabilities — that is, to the 
case T]^-^- <^ 1. In this case, we can drop the exponential and the probabilities of observing 
one photon at c or d, given the observation of one photon at c are: 

P,,0|,0,^,(^^.-^), (32) 

.,oaM)^,(^-^). ,33, 

In particular, if a = 6, then Eq. ()32|) and ()33p imply that, after the detection of a photon at 
c, the probability of subsequently detecting another c photon is 3 times the probability of 

9 



detecting a d photon. This, of course, requires that the times between detections be short 
compared to the period corresponding to the frequency difference of the two lasers. 



2. Multiple Detections 

Muhiple detections can determine the phase completely. Using a naive expression like 
Eq. (fTT|). one could argue that after Mc detections at c and Ma detections at d, in the limit 
of large M = + M^, the probabilities of detecting Nc photons at c or photons at d 
are: 

hm F(iVe|M„ M,) = l-i + a6cos0o e-^(^+'^^™^*») (34) 



and 



hm P(iV,|M„M,) = |^f^^-a6cos0o) ^ (35) 

where cos 0o is given by the expression 

f a^ + b'M,-M, \ 
'''^^=[ 2ab M^ + mJ- 

Note that Eq. ()36|) determines <po only up to a sign, but — since Eqs. ()34p and (jH^I) only 
depend on cos (po — the results are unambiguous. The result of Eq. (jHBj) is only obtained if 
the detection efficiency is the same at c and d. The probabilities of Eqs. fl34|) and ()35|) are 
the same results that would follow from coherent states with a fixed phase difference 0o- 
Effectively, the multiple detections have determined the phase (difference) completely. This 
is an illustration of what Ref. calls a "phase-lock without phase" . 

One could worry that Eq. (jHUj) gives a complex phase if 
a-priori probability of detecting and photons is 



Mc-Ma 



Mc+Ma 



> -J^. However, the 



P(M„M,) = ^^y^ ^(^— + a6cos0)^^^^(^— -a6cos0)^"<'e-''^'^+''^ (37) 

Now, the function to be integrated in Eq. has, for large Mc and M^, sharp peaks at 
cos 00 given in Eq. (jHUj) and is neglegible elsewhere. So the probability of Eq. (jHTj) will vanish 



Mc-Mg 

Mc+Md 



< 



^2 ° j,2 , because otherwise the peaks would be outside the 



for large M^, unless 
integration region. 

It should be obvious that it is possible to rigorously derive the results of Eqs. (jH^ and 
(jnH) only if Mc + M^ ^ Nc, N^. Otherwise, Eq. represent the only prediction for the 
future detections that follows from the initial observation. This means that it is not entirely 
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correct to state as in Ref. js^ that "appropriate measurements of part of a laser beam will 
reduce the quantum state of the rest of the laser beam to a pure coherent state" . The phase 
measurement on the first package will "reduce" only the first package: the phase of the rest 
of the beams will be "determined" only in a statistical way. Again, this is a "weakness" of 
the Bayesian methods used. In practice, however, a statistical state reconstruction could be 



sufficient 



of Ref. [5] are still correct in practice. 
B. Different Frequencies 

I will now briefly discuss the case in which the two lasers have different frequencies. Up 
to this point, we have assumed that the time length of the photon detections and the time 
separation of initial and flnal detections are short with respect to the period corresponding 
to the difference of the frequencies of the two lasers. If only the length of the detection 
processes is short, but the time difference between the photon detections is not, then the 
probability of detecting photons at c, at time t — after detecting Mc photons at c and 
photons at d, at time — is 

P{N,\M,,M,) = ^ / rf0P(0|M„M,)(^— + a6cos(0-Aa;t))^=e-''(^+'^''^°'^(*-^-*«, 

^^c- Jo ^ 

(38) 

where Au is the frequency difference: 

Au = Ua - LJb, (39) 

and P{(j)\Mc, Md) is given by Eq. (jHUI) (Eq. is derived by substituting the value 

»c^at^j)e"^fe '"fe* g^^^^ — substitutions made possible by the normal ordering.) 
Because of the ambiguity in the sign of in Eq. (j3Up . measurements at different times are 
necessary to get a complete picture. If one earlier detects Mdti) photons at c, and Mditi) 
photons at d at time ti — where i = 0,1,. ..m — the conditional probability of detecting Nc 
photons at c, at time t is 

P{N,\M,{U),Md{U)) = 



y accurate for practical applications, like quantum cryptography, so the conclusions 



4- / d0P(0|M,(tO,M,(t,))(^— + a6cos(0-Aa;t))^^e-^(^+'^''=°^(*-^-*», (40) 
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were P{4>\Mc(ti), Md{ti)), the conditional probability of the phase having the value (j), 
after the measurements at times to,ti . . .tm, is given by the expression 



P{<t>\M,{U),Md{U)) = 

UT=o [i^ + abcos{cl)-Auti)y'^^''\^-abcos{(j)-Auti))^'^^'^ 



(41) 



JT ^^0' niio ( ( ^ + cos (0' - Au;ti ) ) ) ( - a6 cos (0' - Acufi ) ) ) ) 

Eqs. fj40p and 1)411) describe how coherent oscillations can be predicted by earlier measure- 
ments. They give an analytic description of the process of "generating" coherent oscillations 
from the entanglement produced by photon detections, as described by Klaus M0lmer in 
Ref. Eqs. (pn|) and (jUj) also provide a quantitative prediction from the model of Ref. 
for propagating laser beams, a prediction that could be used to test the model. 

V. CONCLUSION 

The discussion presented in the previous section is, hopefully, a contribution to the con- 
tinuing discussion of the problem of quantum phase. Differently from other authors (see, 
for instance, Ref. 0]), I did not try to define a phase "observable". In my approach, phase 
is a parameter of an (initially) unknown quantum state. Photon detections determine a 
probability distribution for the unknown phase. 

The information-theoretic approach to quantum phase used in this paper has a precedent 



in Ref. jlOj] (which, however, dealt with neutron interferometry) . The approach used in 
the present paper does not make use of semiclassical methods and takes into account both 
photon statistics and the finite efficiency of the detectors. 
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